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High-precision approximate analytic expressions for energies and wave functions are found for 
arbitrary physical potentials. The Schrodinger equation is cast into nonlinear Riccati equation, 
which is solved analytically in first iteration of the quasi-linearization method (QLM). The zeroth 
iteration is based on general features of the exact solution near the boundaries. The approach is 
illustrated on the Yukawa potential. The results enable accurate analytical estimates of effects of 
parameter variations on physical systems. 



We find accurate analytic presentation of wave func- 
tions and energies for an arbitrary physical potential 
U{r). We use the quasilinearization method (QLM) 
suggested recently for solving the Schrodinger equa- 
tion after conversion to Riccati form In QLM 
the nonlinear terms of the differential equation are ap- 
proximated by a sequence of linear expressions. The 
QLM is iterative but not perturbative and gives stable 
solutions to nonlinear problems without depending on 
the existence of a smallness parameter. 



Substitution of expression y{r) 
radial Schrodinger equation 
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converts the 
U{r)]x{r) = 



Ex(r) into the nonlinear Riccati equation y'{r) + 
y\r) =k^ir),k'^{r) = 2m[U{r)-E]. 

The corresponding QLM equation is y'n+i{r) + 
2j/„+iy„(r) = yl{r) + kl{r), where kl is obtained 
from fc^(r) by replacing there E by energy of n- 

/ Xi^ (^)-f^{^)XTi (^)f^^ 

th iteration E^ = ■'° . . — — — ,Xn{r) = 

Jg X„{nXn{r)dr 

C exp (a yn{r')dr'^ . Since the_QLM iterations have 
very fast quadratic convergence 1, 2], one can expect 
that even the first iteration which is given by an ana- 
lytic expression 
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will be accurate if the zeroth iteration based on gen- 
eral features of solutions near the boundaries is cho- 
sen. 

For illustration we find the wave functions and bind- 
ing energies of Yukawa potential analytically in the 
first QLM iteration. To estimate the precision of our 
analytic solution we solve the Schrodinger equation 
numerically as well. 

The Yukawa potential U{r) = —g- , g > Q 

was suggested in the early days of quantum mechan- 
ics for description of nucleon interactions. During 
last decades the Yukawa potential have been used 
in atomic physics applications, such as the screening 
of nucleon electromagnetic field by electron cloud, or 
atoms under external pressure and in connection to 
quark interactions with parameters A and g depend- 
ing on the temperature of the quark-gluon plasma. 



Let us try to guess the simplest form of the wave 
function. The large distance behavior is ipir) ^ e~^^ 
with rj — \/ —2mE, while the small distance behavior 
is determined by the Kato condition Q "0(0) = ~ ^ 
with /i = mg. Noting also that the radial wave func- 
tion should have a nonzero value at the origin, we 
come to the following initial guess function V'o('') ~ 

^ with a chosen to satisfy the Kato condition 

which leads to a = 2/i — r;. Thus we find for the initial 
guess Xo('') = W-'o('') = [e"'''' — e"^^''"'')''] , and 
therefore yaif) = — /i + (/i — yy) coth [{fi — rj) r] where 

is the normalization factor. 



Inserting this into the equation for Eq and using a 
straightforward integration, one obtains for the zeroth 
order ground state energy 
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Since Eq{t^) = — ^ this is a transcendental equa- 
tion for the parameter 7/. 

The first iteration of the logarithmic derivative is 



J 



given by Eq. ([T]). Its explicit form is given by yi{r) 
Voir) + where (3,] 
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TABLE I: The energy values calculated by direct numeri- 
cal solution of the wave equation (Ed), and in the zeroth 
and first iterations of QLM (-Eo) and (-Ei). 



A 


-Eq 


-El 


-Ed 


0.2 


0.32679 


0.32680851 


0.32680851 


0.5 


0.14795 


0.1481170 


0.1481170 


0.8 


0.04445 


0.0447042 


0.0447043 



Here Ei{z) = — ^^dt is exponential integral func- 
tion. Inserting yi{r) into the expression for Ei, one 
can calculate the first iteration energy. 

Numerical results for the binding energies are given 
in Table 1 while the typical modified ground state 



wave function x(r) is displayed on Fig. 1. We are 
using atomic system of units n2 ^ g ~ I. The di- 
mensionless parameter A is expressed in units of the 
inverse Bohr radius while the energy E is expressed 
in Hartree. 

Comparison of our approximate analytic expres- 
sions for binding energies and wave functions with the 
exact numerical solutions demonstrates their high ac- 
curacy in the wide range of physical parameters. The 
accuracy ranging between 10"*^ and 10~^ for the ener- 
gies and, correspondingly, lO^'^ and 10~* for the wave 
functions is reached. Similarly, the accurate analytic 
presentation of the solution of the Schrodinger Equa- 
tion could be obtained for an arbitrary potential. 
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FIG. 1: Ground state wave functions for the Yukawa po- 
tential with parameter A = 0.2. The exact modified wave 
function (sohd fine) is depicted at the upper part of the 
graph. Relative logarithmic deviations from the exact val- 
ues for the QLM wave functions of the zeroth (dash line) 
and the first order (dot line) are shown in the lower part. 
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